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[57] ABSTRACT 

Disclosed is a multiple measurement multiple energy 
X-ray imaging system in which a plurality of measure- 
ments are processed to provide a first image signal rep- 
resenting a desired parameter of an object and in which 
the plurality of measurements are processed to provide 
a second processed image signal having greater signal- 
to-noise ratio than the first processed image signal but in 
which extraneous artifacts may be introduced into the 
signal. The spatial location of edges of the extraneous 
artifacts are determined. The first processed image sig- 
nal and the second processed image signal are combined 
to provide an improved image signal except at the spa- 
tial locations of the extraneous artifacts where the first 
processed image signal is used for the image signal 
thereby eliminating the extraneous artifact from the 
displayed image. 

14 Claims, 4 Drawing Figures 
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improved selective image is formed by combining the 
MULTIPLE MEASUREMENT NOISE REDUCING selective image with the lower-noise image using the 
SYSTEM USING ARTIFACT EDGE lower frequency components of the former and the 

IDENTIFICATION AND SELECTIVE SIGNAL higher frequency components of the latter. To assure 
PROCESSING 5 the proper ampltitude of the higher frequency compo- 

nents, the components are weighted such as with the 

This invention relates generally to imaging systems ratio of the derivative of the selective image to the 
such as X-ray systems, and more particularly the inven- derivative of the lower-noise image. Nonlinear thresh- 
tion relates to the processing of images derived from a olds are used for those cases where this ratio becomes 
plurality of measurements. io an unreliable indicator of the correct amplitude of the 

In many imaging application areas, images are con- high frequency components. One limitation of the 
structed as a weighted sum of a plurality of measure- Macovski system is the introduction of extraneous sig- 
ments. A prime example is the recent new develop- nals (herein termed "artifacts*') by the low noise signal, 
ments in X-ray imaging. Here measurements are made The present invention is directed to an improvement 
at different energies and/or at different times. These 15 in the Macovski system. The improvement eliminates 
measurements are then combined to provide selective severe artifact edges which may not otherwise be sup- 
images representing specific materials. pressed in the Macovski system. Briefly, the locations of 

One example is that of temporal subtraction tech- the artifact edges in the image signal are identified and 
niques using digital radiography as described in the the display signal is switched from the Macovski noise- 
publication by C. A. Mistretta and A. B. Crummy, "Di- 20 reducing system (NRS) signal to the initial selective 
agnosis of Cardiovascular Disease by Digital Subtrac- signal, S, or a signal derived therefrom. In a preferred 
tion Angiography," in Science, Vol. 214, pp. 761-65, embodiment, the artifact edges are identified by passing 
1981. Here measurements taken before and after the an artifact image signal through a plurality of derivative 
administration of iodine into blood vessels are sub- operators in the horizontal and vertical directions to 
tracted to provide an image of vessels alone. Another 25 provide the two components of the spatial gradient 
example is that of energy-selective radiography as de- vector, then taking the magnitude of the gradient at 
scribed in the publication by L. A. Lehmann, et al., each point. This magnitude is then applied to a thresh- 
"Generalized Image Combinations in Dual KVP Digi- old detector. Since severe edges are of interest, a rela- 
tal Radiography," in Medical Physics, Vol. 8, pp. tively large threshold is used. When the threshold is 
659-67, 1981. Here measurements made at different 30 exceeded, thereby indicating a severe artifact edge in 
energies are combined to enhance or eliminate specific the NRS signal, the selective signal, S, or a signal de- 
materials. A third example is energy-selective comput- rived therefrom is applied to the display, 
erized tomography as described in the publication by R. Accordingly, an object of the present invention is a 
E. Alvarez and A. Macovski, "Energy-Selective Re- multiple measurement noise reducing system in which 
constructions in X-ray Computerized Tomography," in 35 severe artifacts are eliminated. 

Physics in Medicine & Biology, Vol. 21, pp. 733-744, Another object of the invention is a method of im- 
1976. Here sets of measurements are made at two en- proving the image in a multiple measurement noise 
ergy spectra and distortion-free reconstructions are reducing system. 

made which delineate the average atomic number and A feature of the invention is the detection of artifact 
density of each region. A fourth example is hybrid sub- 40 edge locations in an image signal and switching to a 
traction described in U.S. application Ser. No. 260,694, signal not having the artifacts at the detected locations, 
filed May 5, 1981 and in the publication by W. R. The invention and objects and features thereof will be 
Brody, "Hybrid Subtraction for Improved Intravenous most readily apparent from the following description 
Arteriography," in Radiology, Vol. 141, pp. 828-831, and appended claims when taken with the drawings, in 
1981. Here dual energy measurements are made before 45 which: 

and after the administration of iodine. Each dual energy FIG. 1 is a block diagram of an X-ray system for 
pair is used to eliminate soft tissue. In this way the resul- acquiring signals useful in the invention, 
tant subtracted vessel images are immune to soft tissue FIG. 2 is a block diagram of an X-ray system in ac- 
motion- cordance with the prior art. 

In each case, where a number of measurements are 50 FIG. 3 is a functional block diagram of an X-ray 
combined to select specific material properties, the re- signal processing system in accordance with the present 
sultant SNR Signal-to-noise ratio) is reduced as com- invention. 

pared to the nonselective image. Therefore, these vari- FIG. 4 is a functional block diagram of edge detec- 
ous processing techniques which improve the visualiza- tion circuitry useful with the system of FIG. 3. 
tion of disease processes by selecting specific materials 55 Referring now to the drawings, FIG. 1 is a block 
have the disadvantage of a reduction in SNR. This diagram of an X-ray system which includes an X-ray 
reduced SNR can interfere with the ability to visualize tube 8 supplied by variable voltage power source 3. The 
regions of disease. X-rays are transmitted through object 1 F normally the 

A system to improve the SNR of the hybrid subtrac- human anatomy. The transmitted rays are received by 
tion and other multiple measurement systems has been 60 X-ray detector 4, such as an image intensifier and televi- 
developed by A. Macovski and is disclosed in U.S. sion camera, or a one- or two-dimensional array of 
patent application Sen No. 415,334 filed Sept 7, 1982 scintillators, etc. The image signal 5 is applied to storage 
for "Multiple Measurement Noise Reducing System." system 6 where a plurality of images 9, 10, 11, 12, etc., 
An array of measurements are taken of an object under can be stored. 

different conditions. These are combined to select a 65 In temporal subtraction a first image 9 is stored. Fol- 
specific aspect of the object. The same measurements lowing adminstration of a contrast agent using syringe 
are then combined in different weightings to provide a 2, a second image 10 is stored. These are then used in 
lower-noise image, without the desired selectivity. The the subsequent processing systems. In energy selective 
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systems, images 9, 10, 11, 12, etc., are stored at different 
values of beam energy corresponding to different anode 
voltages 3. In dual-energy systems, two voltages are 
used. In addition different X-ray beam filters, not 
shown, can be added. 5 

In the Brody hybrid subtraction system previously 
described, measurements at two voltages are taken prior 
to the adminstration of contrast material. These are 
stored as 9 and 10. Following administration of the 
contrast agent, two additional measurements at the 10 
same two voltages are taken and stored as 11 and 12. 
These four measurements are then processed to obtain a 
low-noise image of the vessels opacified by the contrast 
material from syringe 2. 

Referring now to FIG. 2, the lines 9, 10, 11 and 12 15 
representing an array of X-ray measurements of an 
object, are processed in accordance with the Macovski 
system described above. These signals could be ob- 
tained, for example, from the output of the television 
camera in a fluoroscopic system where an array of im- 
ages are formed and stored, for example, on a magnetic 
disc or in a solid-state frame memory. These signals can 
represent X-ray images taken with different photon 
energy spectra or taken before and after the administra- 
tion of a contrast agent. 

In general the multiple measurements are taken so as 
to provide the isolation or elimination of specific mate- 
rials. One simple example involves the isolation of ves- 
sel images by subtracting measurements taken before 3Q 
and after the administration of an iodinated contrast 
agent. Other examples include measurements at differ- 
ent X-ray photon energy spectra which are then com- 
bined to eliminate or enhance specific materials. 

Weighted summer 13 takes the measurements and 35 
provides selective signal 17, containing the desired ma- 
terial properties. Unfortunately, however, signal 17 
often has relatively poor SNR. The weighted sum usu- 
ally involves negative weights which diminish the sig- 
nal. The noise variances of each independent measure- 43 
ment add, however, to provide a large variance, hence 
a low SNR. This SNR is considerably improved using 
low pass filter 15 which reduces the bandwidth. This 
noise reduction, however, is accomplished at the price 
of reduced resolution. A primary objective of this in- 45 
vention is to restore the resolution using a source of 
high frequency information with reduced noise. 

Signals 9, 10, 11 and 12 are also applied to weighted 
summing structure 14. These weights can be chosen to 
maximize the signal-to-noise ratio for a particular com- 50 
ponent of the image. For example, for vessel imaging 
they can be chosen to maximize the signal-to-noise ratio 
for the iodinated contrast agent. This is in sharp distinc- 
tion to weighted summing structure 13 where the 
weights are chosen to isolate iodine and eliminate the 55 
various tissue components. Ideally, the weights which 
maximize the SNR will vary in different portions of the 
image, depending on the intervening tissue. Thus, the 
weights can be varied as the signals are scanned. How- 
ever, for simplicity, the weights can be fixed to maxi- 60 
mize the average SNR with relatively small degrada- 
tion in performance. 

The high SNR weighted-sum signal, 18, is applied to 
high pass filter 16 which extracts the high frequency 
components. A preferred embodiment of this high pass 65 
filter is the complement of low pass filter 15. Thus the 
sum of the normalized transmissions equals unity within 
the frequencies of interest. These filters can either be 



implemented in the frequency domain or as convolution 
functions. 

Combiner 21 combines the low frequency image or 
signal 19 with the high frequency image or signal 20 to 
obtain the processed image 22 which is displayed in 
display 23. Since signal 19 is low noise because of the 
filtering, and signal 20 is low noise because of the 
weightings, the processed signal 22 is a low noise ver- 
sion of the desired signal, having the full bandwidth or 
resolution. 

The degree of artifacts in the image can depend on 
the nature of the combiner 21. For simplicity, the com- 
biner 21 can simply be an adder. Here the desired se- 
lected image will have the required low noise and high 
resolution. However, signal 20 contains the high fre- 
quency components of undesired structures. For exam- 
ple, in vessel images, signal 20 contains the edges of 
undesired structures such as bone and various soft tissue 
organs and motion artifacts. Since these are only edges, 
however, a trained radiologist can ignore them and 
concentrate on the vessels; 

This system is particularly applicable to the previous- 
ly-described hybrid subtraction scheme. Assume sum- 
mer 13 weights the four signals to eliminate everything 
but iodine, while also subtracting soft tissue motion 
artifacts. Summer 14 can, for example, weight the sig- 
nals to provide a low noise temporal subtraction, with- 
out concern about soft tissue motion artifacts. This can 
be accomplished by using a large positive weight for the 
low energy signal before iodine is administered, and a 
smaller positive weight for the high energy signals since 
it contains less of the iodine component. These weights 
are reversed for the measurements taken after the iodine 
administration to provide a temporal subtraction of 
everything except iodine. 

In this case, in the absence of motion, when 19 and 20 
are added, the desired isolated vessel signal will be 
formed having its full resolution. In the presence of soft 
tissue motion, signal 19 will be immune and continue to 
be vessels only, as described in the previously refer- 
enced application on hybrid subtraction. Low noise 
signal 20, however, will contain motion artifacts. Thus 
when they are added in combiner 21, the edges of these 
artifacts will be present. 

In accordance with the present invention, the multi- 
ple measurement noise reducing system of Macovski is 
modified whereby severe artifact edges, which may not 
be reliably prevented in the Macovski system, are elimi- 
nated by switching to the parameter selective signal or 
a signal derived therefrom in the image regions where 
the severe artifact edge is present. 

Apparatus in accordance with the invention is illus- 
trated in the block diagram of FIG. 3. The selective 
signal, S. 17 and the low noise signal, M, 18 ar applied 
to the NRS filter and combiner shown generally at 30 
and which correspond to the functional elements 15, 16. 
and 21 in FIG. 2. The low noise signal, M, or a signal, 
D, which is the difference (M — S) or any other signal 
which includes the artifact, is applied to an edge detec- 
tion operator 32 which identifies artifact edges and 
generates a signal, £, in response thereto. Tbe signal E 
is applied to a switch 36 whereby the signal from the 
NRS filter and combiner 30 is overridden upon detec- 
tion of an artifact edge, and the selective signal. S. or a 
spatially filtered version thereof derived from the spa- 
tial filter 34, is applied as the output image signal, S. 
Accordingly, in the region of the artifact edge, as identi- 
fied by the edge detection operator, the artifact gener- 
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ated in the NRS filter and combiner 30 is eliminated by 
disconnecting the NRS filter and combiner. 

FIG. 4 is a functional block diagram of one embodi- 
ment of an edge detector operator 32 of FIG. 3 and 
includes a plurality of derivative operators (or other 
. edge filters) 40, 41, and 42, each receiving the input 
image signal, D or M, in the horizontal and vertical 
directions to provide the two components of the spatial 
gradient vector. The magnitude of the gradient vector 
at each point in the image signal as derived by the point 
operator 44 is then compared with a threshold. Other 
edge map techniques can be employed such as template 
filters known in the art to detect edge orientation and 
strength. 

There has been described an improved multiple mea- 
surement noise reducing system in which the presence 
of severe artifacts is eliminated by detecting the pres- 
ence of artifacts and selectively overriding the noise 
reduction system (NRS). By detecting the edges of an 
artifact through use of spatial operators and then select- 
ing between the NRS signal and a higher noise or 
slightly distorted signal without artifacts, an improved 
image is provided. While the invention has been de- 
scribed with reference to a specific embodiment, the 
description is illustrative of the invention and is not to 
be construed as limiting the invention, Various modifi- 
cations and applications may occur to those skilled in 
the art without departing from the true spirit and scope 
of the invention as defined by the appended claims. 

What is claimed is: 

1. In an imaging system, apparatus for reducing extra- 
neous image signals in a processed image signal of an 
object representing a desired parameter derived from a 
plurality of measurements comprising: 

first processing means for processing said plurality of 
measurements and producing a first processed 
image signal representing said desired parameter, 

second processing means for processing said plurality 
of measurements and producing a second pro- 
cessed image signal having a greater signal-to-noise 
ratio from said first processed image signal, 

third processing means interconnected with said first 
processing means and said second processing 
means for receiving and combining said first pro- 
cessed image signal and said second processed 
image signal and producing a third processed 
image signal, 

means for determining the spatial location of edges of 
extraneous images in said second processed image 
signal, and 

means for switching from said third processed image 
signal to a signal based on said first processed 
image signal when edges of extraneous images are 
detected. 

2. Apparatus as defined by claim 1 wherein said 
means for determining the spatial location of edges 
includes differential means for determining the compo- 
nents of the spatial gradient signal vector of the differ- 
ence between the said first processed image signal and 
said second processed image signal. 

3. In an imaging system, apparatus for reducing extra- 
neous image signals in a processed image signal of an 
object representing a desired parameter derived from a 
plurality of measurements comprising: 

first procssing means for processing said plurality of 
measurements and producing a first processed 
image signal representing said desired parameter, 



second processing means for processing said plurality 
of measurements and producing a second pro- 
cessed image signal having a greater signal-to-noise 
ratio from said first processed image signal, 
first filter means for receiving and low pass filtering 

said first processed image signal, 
second filter means for receiving and high pass filter- 
ing said second processed image signal, 
means for combining said filtered first processed 
image signal and said filtered second processed 
image signal to produce a combined image signal, 
means for determining the spatial location of edges of 
extraneous images in said second processed image 
signal, and 

means for switching from said combined image signal 
to a signal based on said first processed image sig- 
nal in response to detection of edges of extraneous 
images. 

4. Apparatus as defined by claim 3 wherein said 
means for determining the spatial location of edges 
includes differential means for determining the compo- 
nents of the spatial gradient signal vector of the differ- 
ence between said first processed image signal and said 
second processed image signal. 

5. Apparatus as defined by . claim 3 wherein said first 
filter means and said second filter means are comple- 
mentary in frequency. 

6. Apparatus as defined by claim 3 wherein said first 
processed image signal is derived from a first weighted 
sum of said measurements and said second processed 
image signal is derived from a second weighted sum of 
said measurements. 

7. Apparatus as defined by claim 6 wherein said imag- 
ing system comprises an X-ray system and said measure* 
ments are at different X-ray energies. 

8. Apparatus as defined by claim 6 wherein said imag- 
ing system comprises an X-ray system and said measure- 
ments are at different times. 

9. For use in an X-ray imaging system in which a 
plurality of measurements are made, a method of reduc- 
ing extraneous signals in a processed image signal com- 
prising the steps of: 

processing said plurality of measurements to provide 
a first processed image signal representing a de- 
sired parameter of an object, 
processing said plurality of measurements to provide 
a second processed image signal having greater 
signal-to-noise ratio than said first processed signal 
and introducing artifact images, 
combining said first and second processed image 
signals and thereby providing a third processed 
image signal, 

determining the spatial location of edges of said arti- 
facts in said second processed signal, and 
switching from said third processed image signal to 
said first processed image signal at the spatial loca- 
tion of said artifacts. 

10. For use in an imaging system in which a plurality 
of measurements are made, a method of reducing extra- 

60 neous signals in a processed image signal of an object 
comprising the steps of: 

processing said plurality of measurements to provide 
a first processed image representing a desired pa- 
rameter of said object, 
processing said plurality of measurements to provide 
a second processed image signal having greater 
signal-to-noise ratio than said first processed signal 
and introducing artifacts, 
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low pass filtering said first processed image signal to 
reduce its noise, 

high pass filtering said second processed image signal, 
said high pass filtering being complementary to 
said low pass filtering, 5 

combining said filtered first processed image and said 
filtered second processed image to provide an im- 
proved processed image signal, 

determining the spatial location of edges of said arti- 
facts, and 10 

selectively switching from said improved processed 
image signal to a signal based on said first pro- 
cessed image signal at the spatial locations of arti- 
facts. 

11. The method as defined by claim 10 wherein said 15 
step of determining the spatial location of extraneous 



S 

artifact edges include determining the components of 
the spatial gradient signal vector of the difference be- 
tween said first processed image signal and said second 
processed image signal. 

12. The method as defined by claim 10 wherein said 
first processed image signal is derived from a first 
weighted sum of said measurements and said second 
processed image signal is derived from a second 
weighted sum of said measurements. 

13. The method as defined by claim 10 wherein said 
imaging system comprises an X-ray system and said 
measurements are at different X-ray energies. 

14. The method as defined by claim 10 wherein said 
imaging system comprises an X-ray system and said 

measurements are at different times. 

• » « * « 
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